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ABSTRACT 

We construct a maximum-likelihood algorithm - MAXLIMA, to derive the mass 
distribution of the extrasolar planets when only the minimum masses are ob- 
served. The algorithm derives the distribution by solving a numerically stable 
set of equations, and does not need any iteration or smoothing. Based on 50 
minimum masses, MAXLIMA yields a distribution which is approximately flat in 
log M, and might rise slightly towards lower masses. The frequency drops off very 
sharply when going to masses higher than 10 Mj, although we suspect there is 
still a higher mass tail that extends up to probably 20 Mj. We estimate that 5% 
of the G stars in the solar neighborhood have planets in the range of 1-10 Mj with 
periods shorter than 1500 days. For comparison we present the mass distribution 
of stellar companions in the range of 100-lOOOMj, which is also approximately 
fiat in log M. The two populations are separated by the "brown-dwarf desert" , 
a fact that strongly supports the idea that these are two distinct populations. 
Accepting this definite separation, we point out the conundrum concerning the 
similarities between the period, eccentricity and even mass distribution of the 
two populations. 

Subject headings: binaries: spectroscopic — methods: statistical — planetary 
systems — stars: low-mass, brown dwarfs — stars: statistics 



1. INTRODUCTION 

Since the detection of the first few extrasolar planets their mass distribution was recog- 
nized to be a key feature of the growing new population. In particular, the potential of the 
high end of the mass distribution to separate between planets on one side and brown dwarfs 
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and stellar companions on the other side was pointed out by numerous studies (e.g., Basri 
& Marcy 1997; Mayor, Queloz & Udry 1998; Mazeh, Goldberg & Latham 1998). A clear 
mass separation between the two populations could even help to clarify one of the very basic 
questions concerning the population of extrasolar planets — the precise definition of a planet 
(Burrows et al. 1997; see a detailed discussion by Mazeh & Zucker 2001). The lower end of 
the mass distribution could indicate how many Saturn- and Neptune-like planets we expect 
on the basis of the present discoveries, a still unsurveyed region of the parameter space of 
extrasolar planets. 

The present number of known extrasolar planets — more than 60 arc known (Encyclope- 
dia of extrasolar planets, Schneider 2001), offers an opportunity to derive a better estimate of 
the mass distribution of this population. In order to use the derived masses of the extrasolar 
planets we have to correct for two effects. The first one is the unknown orbital inclination, 
which renders the derived masses only minimum masses. The second effect is due to the 
fact that stars with too small radial-velocity amplitudes could not have been detected as 
radial-velocity variables. Therefore, planets with masses too small, orbital periods too large, 
or inclination angles too small are not detected. 

The effect of the unknown inclination of spectroscopic binaries was studied by numerous 
papers (e.g., Mazeh & Goldberg 1992; Hcacox 1995; Goldberg 2000), assuming random 
orientation in space. Heacox (1995) calculates first the minimum mass distribution and then 
uses its relation to the mass distribution to derive the latter. This calculation amplifies 
the noise in the observed data, and necessitates the use of smoothing to the observed data. 
Mazeh & Goldberg (1992) introduce an iterative algorithm whose solution depends on the 
initial guess. In the present work we followed Tokovinin (1991, 1992) and constructed an 
algorithm — MAXimum Likelihood MAss, to derive the mass distribution of the extrasolar 
planets with a maximum likelihood approach. MAXLIMA assumes that the planes of motion 
of the planets are randomly oriented in space and derives the mass distribution directly by 
solving a set of numerically stable linear equations. It does not require any smoothing of 
the data nor any iterative algorithm. MAXLIMA also offers a natural way to correct for the 
undetected planets. 

The randomness of the orbital planes of the discovered planets were questioned recently 
by Han et al. (2001), based on the analysis of Hipparcos data. However a few very recent 
studies (Pourbaix 2001; Pourbaix & Arenou 2001; Zucker & Mazeh 2001a,b) showed that 
the Hipparocos data do not prove the nonrandomness of the orbital planes, allowing us to 
apply MAXLIMA to the sample of known minimum masses of the planet candidates. 

In the course of preparing this paper for publication we have learned about a similar 
paper by Jorissen, Mayor & Udry (2001) that was posted on the Astrophysics e-Print Archive 
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(astro-ph). Like Heacox (1995), Jorrisen et al. derive first the distribution of the minimum 
masses and then apply two alternative algorithms to invert it to the distribution of planet 
masses. One algorithm is a formal solution of an Abel integral equation and the other is 
the Richardson-Lucy algorithm (e.g., Heacox 1995). The first algorithm necessitates some 
degree of data smoothing and the second one requires a series of iterations. The results of 
the first algorithm depend on the degree of smoothing applied, and those of the second one 
on the number of iterations performed. MAXLIMA has no built-in free parameter, except the 
widths of the histogram bins. In addition, Jorissen et al. did not apply any correction to the 
selection effect we consider here, and displayed their results on a linear mass scale. We feel 
that a logarithmic scale can illuminate some other aspects of the distribution. Despite all 
the differences, our results arc completely consistent with those of Jorrisen et al., the sharp 
cutoff in the planet mass distribution at about 10 Mj, and the small high-mass tail that 
extends up to about 20 Mj in particular. 

Section 2 presents MAXLIMA, while section 3 presents our results. Section 4 discusses 
briefly our findings. 



Our goal is to estimate the probability density function (PDF) of the secondary mass 
— /m(m)cim, given a set of observed minimum masses {yj',j = 1, A^}, where yj = rrij x sinij, 
rrij is the mass of the j-th planet and ij is its inclination. Within certain assumptions and 
limitations, MAXLIMA finds the function fm that maximizes the likelihood of observing these 
minimum masses. 

Note that the present realization of MAXLIMA uses the approximation that the mass of 
the unseen companion is much smaller than the primary mass. Within this approximation 
yj can be derived from the observations of each system, given the primary mass. In general, 
when the secondary mass is not so small, the value of rrij x sin ij cannot be derived from the 
observations, and a more complicated expression has to be used. Nevertheless, the extension 
of MAXLIMA to those cases is straightforward, and will be worked out in details in a separate 
paper. 

We assume that the directions of the angular momenta of the systems are distributed 
isotropically, which will cause sin i to have a PDF of the form: 



2. 



MAXLIMA 



2.1. 



The Unknown Inclination 
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where we denote sin i by s. We further assume that the planet mass and its orbital inclination 
are uncorrelated, and therefore the joint PDF of sin i and the planet mass have the form: 

fms(rn, s) dmds = fm{m) ^ dmds . 

V 1 — s 

Now, define a variable y — m x sini — m x s, which has the PDF 

/ /„w^dm*= / Mp) (^2) . 

J TTT- J a/1 — [y/my 

A sketch of contours of constant y in the m — s parameter space is plotted in Figure 1. 

We wish to estimate fm{'<Ti) in the form of a histogram with K bins, between the limits 
"^min <Tn< iTimax- We thus cousidcr a partition of the interval [m^j„, mmax)'- 

rrimin = mi,m2, • • • , rriK+i = mmax , 

for which the PDF gets the form: 

fm{m) = bk, ii mk<m< ruk+i, for k = 1, . . . , K. 

Note, that the function fmijn) is supposed to be a probability density function, and 
therefore its integral must equal one: 

K 

J]6fcAmfc = l, (1) 

k=l 

where Am^ = rrik+i — rrik for /c = 1, . . . , K. 
The PDF of y then gets the form: 

= y^^fe / — /, i ^drndy , 

where the integrals do not depend on fm{Tn) at all, but only on the intervals borders and y. 

Now we can solve our problem in a maximum- likelihood fashion by finding the set of hkS 
that maximizes the likelihood of the actually observed values - yj. The likelihood function 
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is: 



where 



N 



C(bi,...,bK;yi,...,yN) = H-^f^^^) 

'•'"'=+1 1 yj/m 



TV K 



TV K 
3=1 k=l 



mfe+i < yj 

A,, = <; ^Tnfe "^'^ ^ < "^^-^ (2) 

rmfc+1 1 yj/m V ■ < TUk. 

and the integrals are easily calculated even analytically. 

In the appendix we present an elegant way to find the h^s that maximize log>C directly, 
without any iterations. 



2.2. Simulation 

In order to check the performance of MAXLIMA and its realization we have performed 
several simulations, some of which are presented in Figure 2. In those simulations we gener- 
ated an artificial sample of planets drawn from populations with different PDFs of the planet 
masses, and inclinations oriented isotropically in space. To make the simulation similar to 
the present work we chose the size of each sample to be 50 planets. We assumed no selection 
effects. We then applied MAXLIMA to the simulated sample, the results of which are plotted 
in Figure 2. 

The three examples of Figure 2 clearly show the power of MAXLIMA. 



2.3. Selection Effect: Undetected Planets 

We assume that the sample is constructed of planets with period, P, between Pmin < 
P < Pmax- We further assume that the search for planets discovered all radial- velocity 
variables with amphtude K larger than -fCmin- We have to correct for planets not detected 
because they induce K smaller than the threshold. To do that we note that the amphtude 



6 



can be written as 




1/3 




2/3 



( 




X(P,Mp,Mi,sini) =204 



sm I m s 



-1 



(3) 



The expression Mp x sini is actually our y. For any given value of y and Mi we can 
derive the maximum possibly detected period — -Pmax-detect; given Xmin- This implies that 
if we know the period distribution and we assume that the period is uncorrclatcd to the 
mass distribution, we can estimate for each of the given y's the fraction of planets with 
long periods that were not detected with the same y. This means that to correct for the 
undetected planets with long periods we have to consider each of the j-th detected systems 
as representing some ctj planets. If -Pmax-detect is smaller than Pmax, then aj is larger than 
unity Otherwise aj is equal to unity 

We then can write a new generalized likelihood as: 



The Appendix shows an easy way used by MAXLIMA to find a maximum to log C 

For simplicity we consider only circular orbits. Eccentricity introduces two effects 
(Mazeh, Latham & Stefanik 1996), the first of which is the dependence of K on the ec- 
centricity e. Equation (3) should include an additional factor of (1 — e^)~^/^, which causes 
K to increase for increasing e. The other factor is the dependence of the detection threshold 
Kuiin on e. Our simplifying assumption about the constancy of Xmin throughout the sample 
breaks down when we consider eccentric orbits. This is so because for eccentric orbits the 
velocity variation tends to concentrate around the periastron passage, and therefore Xmin 
increases for increasing eccentricity. These two effects tend to cancel each other (Fischer & 
Marcy 1992), the net effect depends on the characteristics of the observational search. By 
running numerical simulations Mazeh et al. (1996) have found that if the detection limit 
depends on the r.m.s. scatter of the observed radial-velocity measurements, the two effects 
cancel each other for any reasonable eccentricity. We therefore chose not to include the 
eccentricity of the planets in our analysis. 



TV 



£(6i,...,6K;z/i,---,|/iv;ai,---,ajv) = 
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3. ANALYSIS AND RESULTS 

To apply MAXLIMA to the current known sample of extrasolar planets we considered all 
known planets and brown dwarfs as of April 2001. We consider only G- or K-star primaries 
and therefore excluded Gls 876 from the sample. 

Obviously, the present sample in not complete. In particular, not all planets with long 
periods and small induced radial- velocity amplitudes were discovered and/or announced. To 
acquire some degree of completeness to our sample we have decided, somewhat arbitrar- 
ily, to exclude planets with periods longer than 1500 days and with radial- velocity ampli- 
tudes smaller than 40 m/s. The values of these two parameters determine the correction of 
MAXLIMA for the selection effect, for which we assumed a period distribution which is flat in 
log P. This choice of parameters also implies that our analysis applies only to planets with 
periods shorter than 1500 days. We further assumed that the primary mass is 1 M© for all 
systems. 

In order to be consistent with the selection effects and the correction we applied, we 
included in our analysis only planets that were discovered by the high precision radial- velocity 
searches. We had to exclude HD 114762 and similar objects that were discovered by other 
searches (e.g., Latham et al. 1989; Mayor et al. 1992; Mazeh et al. 1996). This does not 
mean that we assume anything about their nature in this stage of the study. Table 1 lists 
all the known planets with G star primaries from the high precision radial-velocity studies. 
Planets excluded from our analysis are marked by an asterisk. We did not take into account 
the known inclination of the planet around HD 209458. All together we are left with 50 
planets. 

The results of MAXLIMA are presented in the lower panel of Figure 3 on a logarithmic 
mass scale. Each bin is 0.3 dex wide, which means about a factor 2 in mass. The value 
of each bin is the estimated number of planets found in the corresponding range of masses 
in the known sample of planets, after correcting for the undetected systems. To estimate 
the uncertainty of each bin we ran 5000 Monte Carlo simulations and found the r.m.s. of 
the derived values of each bin. Therefore, the errors plotted in the figure represent only 
the statistical noise of the sample. Obviously, any deviation from the assumptions of our 
model for the selection effect induces further errors into the histogram, the assumed period 
distribution in particular. This is specially true for the first bin, where the actual number 
of systems is small and the correction factor large. The value of the first bin is sensitive, for 
example, to the assumed Xmin- Assuming Xmin of 50 m/s increased the value of the first bin 
by more than 50%. 

To compare the mass distribution of the planets with that of the stellar secondaries we 
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plot the latter on the same scale in an adjacent panel of Figure 3. We plot here only two bins, 

with masses between 100 and 1000 Mj. We follow the work of Mazeh (1999b) and Mazeh 
and Zucker (2001), and used for those bins subsamples of binaries found by the CfA radial- 
velocity search for spectroscopic binaries (Latham 1985) in the Carney & Latham (1987) 
sample of the high-proper-motion stars (Latham et al. 2000; Goldberg et al. 2000). For the 
smaller bin we used a subsample that included only the Galaxy disk stars (Goldberg 2000), 
and for the larger-mass bin a subsample of this sample that included only primaries with 
masses higher than 0.7 Mq. The values of those two bins were derived with the algorithm 
of Mazeh & Goldberg (1992). 

Note that the upper panel does not have any estimate of the values of the bins with 
masses smaller than 100 Mj. This is so because the CfA search does not have the sensi- 
tivity to detect secondaries in that range. On the other hand, the lower panel does include 
information on the bins below 100 Mj. This panel presents the results of the high- precision 
radial-velocity searches, and these searches could easily detect stars with secondaries in the 
range of, say, 20-100 Mj. We assume that these binaries were not excluded from the vari- 
ous radial-velocity searches at the first place, and further assume that all, or at least most, 
findings of the various research groups corresponding to this range of masses were already 
published. If these two assumptions are correct, then the lower panel does represent the fre- 
quency of secondaries in the mass range of 20-100 Mj. This panel shows that the frequency 
of secondaries in this range of masses is close to zero. The present analysis is not able to tell 
whether this "brown-dwarf desert" extends up to 60, 80 or 100 Mj. 

The relative scaling of the planets and the stellar companions is not well known. The 
spectroscopic binaries come from well defined samples - 577 stars for the lower-mass bin, 
and 312 stars for the higher-mass stars (Goldberg 2000). However, this is not the case for 
the detected planets, specially because the sample of published planets is not complete and 
also because the search samples of the different groups are not well documented in the public 
domain. We assumed, somewhat arbitrarily, that they come from a sample of 1000 stars, 
and scaled the stellar bins accordingly. We also rescaled the stellar bins to account for the 
fact that their bins are larger, and their period range extends up to 3000 days, assuming a 
fiat distribution in log P. Therefore the values of the stellar bins are our best estimate for 
the number of binaries for 1000 stars within a mass range of 0.3 dex, and up to a period of 
1500 days. 

Obviously the relative scaling of the two panels has a large uncertainty. This scaling 
uncertainty is not refiected in the error bars of the higher panel. Nevertheless we think that 
the comparison is illuminating, as will be discussed in the next section. 
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4. DISCUSSION 

The grand picture that is emerging from Figure 3 strongly indicates that we have here 
two distinct populations. The two populations are separated by a "gap" of about one decade 
of masses, in the range between 10 and 100 Mj. Such a gap was already noticed by many 
early studies (Basri & Marcy 1997; Mayor, Queloz & Udry 1998; Mayor, Udry & Queloz 
1998; Marcy & Butler 1998). Those early papers binned the mass distribution hnearly. Here 
we follow our previous work (Mazeh et al. 1998) and use a logarithmic scale to study the mass 
distribution, because of the large range of masses, 0.5-1000 Mj, involved. The logarithmic 
scale has also been used by Tokovinin (1992) to study the secondary mass distribution in 
spectroscopic binaries, and was suggested by Black (1998) to study the mass distribution 
of the planetary-mass companions (see also Mazeh 1999a,b; Mazeh & Zucker 2001; Mayor 
et al. 2001). The gap or the "brown-dwarf desert" are consistent also with the finding 
of Halbwachs et al. (2000), who used Hipparcos data and found that many of the known 
brown-dwarf candidates are actually stellar companions. 

We will assume that the two populations are the planets, at the low-mass side of Figure 3, 
and the stellar companions at the high-mass end of the figure. Interestingly, the mass 
distribution of single stars extends far below 100 Mj (e.g., Zapatero Osorio et al. 2000; 
Lucas & Roche 2000), indicating that the gap separating the two populations of companions 
apparently does not exist in the population of single stars/brown dwarfs. This difference 
probably indicates different formation processes for single and secondary objects. 

The distribution we derived in Figure 3 suggests that the planet mass distribution is 
almost flat in logM over five bins — from 0.3 to 10 Mj. Actually, the figure suggests a 
slight rise of the distribution towards smaller masses. The distinction between these two 
distributions is not possible at this point, when our knowledge about planets with sub- 
Jupiter masses is very limited. At the high-mass end of the planet distribution the mass 
distribution dramatically drops ofi' at 10 Mj, with a small high-end tail in the next bin. 
Although the results are still consistent with zero, we feel that the small vahie beyond 10 
Mj might be real. The dramatic drop at 10 Mj and the small high-mass tail agree with the 
findings of Jorissen et al. (2001). 

Examination of the two panels of Figure 3 suggests that per equal dex range of masses 
the frequency of stellar secondaries is higher than that of the planets by a factor of about 
2. As we emphasized in the previous section, this is a very preliminary result that should 
be checked by future observations. Nevertheless, the frequency of planets is impressive by 
itself. Our results indicate that about 5% of the stars have planets with masses between 1 
and 10 Mj. This is so because the number of multiple planets in this sample is small, so 
the number of planets considered in the figure is about the number of stars found to have 
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one or more planets. If this frequency extends further down the mass axis to Earth masses, 
wc might find that more than 10% of the stars have planets with periods shorter than 1500 
days. 

The analysis presented here raises the question what mechanism can produce flat or 
approximately flat mass distribution of planets up to 10 Mj. What determines the mass 
of the forming planet? The present paradigm assumes that planets were formed out of a 
protoplanetary disk. Is it the mass, density, angular momentum or the viscosity of the disk 
that determined the planet mass? If planets were formed by accreting gas onto a rocky core, 
is the planet mass determined also by the location or evolutionary phase of the formation 
of the rocky core? Any detailed model of planet formation should account for this mass 
distribution. 

The clear distinction between the two populations suggests that planets and stellar 
companions were formed by two different processes. This is so despite the striking similarity 
between the distributions of the eccentricities and periods of the two populations (Heacox 
1999; Stepinski & Black 2000, 2001; Mayor & Udry 2000; Mazeh & Zucker 2001). Even the 
mass distributions of the two populations might be very similar — approximately flat in 
log M. This is still a conundrum that any formation model for planets as well as for binaries 
needs to solve. 

We are indebted to Yoav Benjamini for illuminating discussions. This work was sup- 
ported by the US-Israel Binational Science Foundation through grant 97-00460 and the Israeli 
Science Foundation (grant no. 40/00) 



We want to flnd the maximum likelihood to observe a given set of observed minimum 
masses {yj]j — Ij-^V}. As usual, it is easier to maximize the logarithm of the likelihood 
function: 



where each Aj^. depends on the corresponding yj through Equation (2). 

The 6fe's are not all independent. They are constrained by Equation (1), and therefore 
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we modify our target function by adding a Lagrange multiplier term: 

log>C = ^ aj X log I ^ Ajkbk I + ^ I XI ^k^^k - 1 I ■ 

j=l \k=l J \k=l J 

The optimization is performed by equating the partial derivatives of this target function 
to zero: 

91og£ aj X Aik , , , , 

= J2 Jk /, +AAm. = 0, forfc = l,...,i^ 



dlogC 



K 



k=l 



The parameter A is eliminated quite easily. We first multiply each of the K equations by 
the corresponding hk and then sum them up to get: 

k=l 3=1 2^1=1 k=l 

Changing the order of summation reduces the first term, after a simple manipulation, to 
simply 

N 

Using the constraint reduces the second term to A, and we finally get: 

Neff + X = 0, 

and we can simply set A = —Neff. The K equations we are now left with are: 

^ 1 

We have a set of K non-hnear equations in K variables - b^s. 

An elegant reduction of the complexity of the problem can be achieved if we set K to 
N, by assigning nik = yk ioT k = 1, . . . ,N. Let us also denote (remember that now K — N): 
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The N equations now look like: 



N 

aj X Ajkhj = NeffAm, 



■k , 



(2) 



which is a system of N linear equations in the N variables hj. We can easily solve for them. 
The problem is even more easily solved when we note that the matrix Aj^ is upper-triangular 
and thus the amount of computations needed for the solution is smaller. Furthermore, 
examination of the integrals involved in the calculation of Ajk (Eq. 2) shows that the matrix 
is very close to being diagonally-dominated, and therefore the set of linear equations is 
numerically stable. Having solved for hj we face again a similar system of linear equations 
in order to solve for bk, coming from the definition of hj: 



In this system the matrix is the transposed matrix of the previous system of linear equations 
if the ctj's are all equal to unity. 

Obviously, we wish to estimate the densities of our original intervals. But these are 
easily calculable from the densities bk calculated above, as simple linear combinations. 
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Fig. 2. — Three simulations of MAXLIMA. In all three panels the dashed line shows the 
input distribution and the histogram the results of MAXLIMA. 
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Fig. 3. — The mass distributions of the planets and the stellar companions 
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Table 1. Substellar companions list 



Name 


Msini 


P 


K 




(Mj) 


(days) 


(m s-l) 


*HD 83443 b 


0.16 


29.83 


14 


*HD 16141 


0.215 


75.82 


11 


*HD 168746 


0.24 


6.41 


28 


*HD 46375 


0.249 


3.02 


35 


HD 83443 c 


0.34 


2.985 


56 


*HD 108147 


0.34 


10.88 


37 


HD 75289 


0.42 


3.51 


54 


51 Peg 


0.47 


4.23 


56 


BD -10°3166 


0.48 


3.487 


61 


*HD 6434 


0.48 


22.09 


37 


HD 187123 


0.48 


3.10 


69 


*Gliese 876 c 


0.56 


30.12 


81 


HD 209458 


0.69 


3.52 


86 


V And b 


0.69 


4.617 


71 


HD 192263 


0.787 


24.36 


68 


HD 38529 


0.81 


14.32 


54 


HD 179949 


0.84 


3.09 


101 


55 One 


0.84 


14.65 


77 


*e Eri 


0.86 


2502.1 


19 


*HD 82943 c 


0.88 


222 


34 


HD 121504 


0.89 


64.6 


45 


HD 130322 


1.02 


10.72 


115 


HD 37124 


1.04 


155.7 


43 


p CrB 


1.1 


39.65 


67 


HD 52265 


1.13 


118.96 


45 


*HD 177830 


1.22 


391.6 


34 


HD 217107 


1.27 


7.13 


140 


HD 210277 


1.28 


437 


41 


*HD 27442 


1.43 


426.5 


34 


16 Cyg B 


1.5 


801 


44 


HD 74156 b 


1.56 


51.6 


108 


HD 134987 


1.58 


259.6 


50 


HD 82943 b 


1.63 


445 


46 


*Gliese 876 b 


1.89 


61.02 


210 


HD 160691 


1.97 


743 


54 


HD 19994 


2.0 


454 


45 


HD 213240 


3.7 


759 


91 


V And c 


2.06 


240.6 


58 


HD 8574 


2.23 


228.8 


76 


HR 810 


2.26 


320.1 


67 


47 UMa 


2.39 


1090 


45 


HD 12661 


2.79 


252.7 


88 


HD 169830 


2.96 


230.4 


83 


*14 Her 


3.3 


1654 


73 


GJ 3021 


3.32 


133.82 


164 


HD 92788 


3.34 


326.7 


100 


HD 80606 


3.41 


111.8 


414 


HD 195019 


3.47 


18.2 


272 
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Table 1 — Continued 



Name 


Msini 


P 


K 




(Mj) 


(days) 


(m s-l) 


T Boo 


3.87 


3.31 


469 


Glicsc 86 


4 


15.78 


380 


V And d 


4.10 


1313 


68 


HD 50554 


4.9 


1279 


95 


HD 190228 


5 


1161 


95 


HD 222582 


5.3 


575.9 


184 


HD 28185 


5.6 


385 


168 


HD 10697 


6.35 


1072.3 


119 


HD 178911 B 


6.46 


71.50 


343 


70 Vir 


6.6 


116.7 


318 


HD 106252 


6.81 


1500 


139 


HD 89744 


7.2 


256 


257 


HD 168443 b 


7.2 


58.12 


473 


*HD 74156 c 


>7.5 


2300 


121 


HD 141937 


9.7 


659 


247 


*HD 114762 


11 


84.03 


600 


HD 202206 


14.7 


259 


554 


*HD 168443 c 


15.1 


1667 


288 


*HD 127506 


36 


2599 


891 



